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Abstract 

Dyson-Schwinger equation (DSE) studies show that the 6-quark mass- 
function is approximately constant, and that this is true to a lesser extent for 
the c-quark. These observations provide the basis for a study of the leptonic 
and semileptonic decays of heavy pseudoscalar mesons using a "heavy-quark" 
limit of the DSEs. When exact, this limit reduces the number of independent 
form factors. Semileptonic decays with light mesons in the final state are also 
accessible because the DSEs provide a description of light-quark propagation 
characteristics and light-meson structure. A description of 5-meson decays is 
straightforward, however, the study of decays involving the D-meson indicates 
that c-quark mass-corrections are quantitatively important. 

Dyson-Schwinger equations provide a nonperturbative, Poincare invariant, con- 
tinuum approach to studying quantum field theories: two familiar examples are the 
gap equation in superconductivity and the Bethe-Salpeter equation describing rela- 
tivistic 2-body bound states. As a system of coupled integral equations, a truncation 
of the DSEs is necessary to obtain a tractable problem. The simplest truncation 
scheme is a weak-coupling expansion, which generates every diagram in perturbation 
theory. Hence, in the intelligent application of DSEs to QCD, there is a tight con- 
straint on the ultraviolet behaviour of the Schwinger functions (dressed-propagators 
and vertices). That is crucial in extrapolating into the infrared, in constructing uni- 
formly valid symmetry-preserving truncations, and in developing phenomenological 
models necessary for anticipating the results of the current generation of hadron 
physics facilities. 

The development of efficacious truncations is not a purely algebraic task, and 
neither is it always obviously systematic. Nevertheless, it has become clear ]IJ that 
truncations which preserve the global symmetries of a theory; for example, chiral 
symmetry in QCD, are relatively easy to define and implement and, while it is 
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more difficult to preserve local gauge symmetries, much progress has been made 
with Abelian theories |2| and more is being learnt about non-Abelian ones. In ad- 
dition, contemporary phenomenological applications now address a wide range of 
observables || , yielding qualitatively robust results and a much-needed intuitive un- 
derstanding of many observables inaccessible in perturbation theory. 

A salient feature of the phenomenological application of DSEs is the significant 
role played by the necessary momentum-dependent modification of gluon and quark 
propagators: they are modified in perturbation theory and this modification persists 
and grows in the nonperturbative domain. For example, in a general covariant gauge 
the dressed-gluon propagator is characterised by a single scalar function, which we 
denote T>(k 2 ). Many studies of the gluon DSE show that T>{k 2 ) is strongly enhanced 
in the infrared; i.e, its behaviour in the vicinity of k 2 = can be represented as a 
distribution Q , while for k 2 > 1-2 GeV 2 the perturbative result is reliable. With 
such behaviour manifest in the quark-quark interaction, dynamical chiral symmetry 
breaking (DCSB) and confinement follow without fine-tuning ||. 

These phenomena can be addressed through the DSE for the dressed-quark prop- 
agator: 

1 i 

S (P)'-=- = " 77 2 \ i r>( 2\ -PVvip 2 ) +0- S (j> 2 ) , (1) 

where is the renormalised dressed-quark self energy, which satisfies 

£(p) = {Z 2 -l)i 1 -p + Z, m C + Z x ^ g 2 D iiV {p - q )^- lfl S(q)T a u (q,p), (2) 

with r"(g;p) the dressed-quark-gluon vertex, mf- the current-quark mass, Z\, Z 2 and 
Z 4 renormalisation constants, and ( the renormalisation point. / := J A d 4 g/(27r) 4 
represents mnemonically a translationally-invariant regularisation of the integral, 
with A the regularisation mass-scale. With the infrared-enhanced interaction intro- 
duced in Ref. H and current-quark masses corresponding to 
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one obtains M the dressed-quark mass function depicted in Fig. [I]. It is clear that 
for light quarks (u, d and s) there are two distinct domains: perturbative and 
nonperturbative. In the perturbative domain the magnitude of M(p 2 ) is governed 
by the the current-quark mass, while for p 2 < 1 GeV 2 the mass-function rises sharply. 
This is the nonperturbative domain where the magnitude of M(p 2 ) is determined 
by the DCSB mechanism; i.e., the enhancement in the dressed-gluon propagator. 

For a given flavour, the ratio Cf := Mf/rn} is a single, quantitative measure 
of the importance of the DCSB mechanism in modifying that quark's propagation 
characteristics. As illustrated in Eq. (f|), 
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Figure 1: M(p 2 ) := B(p 2 )/A(p 2 ) obtained in solving the quark DSE. The solution 
of M 2 (p 2 ) = p 2 defines M E , the Euclidean constituent-quark mass. 

this ratio provides a natural classification of quarks as either light or heavy. For 
light-quarks Cf is characteristically 10-100 while for heavy-quarks it is only 1-2. The 
values of Cf signal the existence of a characteristic DCSB mass-scale: M x . At p 2 > 
the propagation characteristics of a flavour with < M x are significantly altered by 
the DCSB mechanism, while for flavours with 3> M x it is irrelevant, and explicit 
chiral symmetry breaking dominates. It is apparent that M x ~ 0.2 GeV~ Aqcd- 

This forms the basis for a simplification of the study of heavy-meson observ- 
ables M that we summarise herein. It motivates an exploration of the fidelity of the 
approximation 

S c ,b(p) = ~ 7— , (5) 



ry • p + M C)b 



where M c ^ ~ 
propagator is 



M c E b , 



so that with p^ := m H v^ 



S c ,b( k + P) 



27 • v 



is the heavy meson velocity, 



ll 

2 k-v-E 
v 2 = -1, 



O 




the heavy-quark 



(6) 



and E > is the difference between 
the heavy-meson mass and the effective- mass of the heavy-quark.) Many simplifi- 
cations follow from neglecting the 1/M-corrections; e.g., it reduces the number of 
independent form factors required to describe heavy-meson — > heavy-meson decays, 
relating them to a minimal number of so-called "universal" form factors, which is a 
characteristic feature of "heavy-quark" symmetry M. The magnitude of M b E makes 
it likely that Eq. (||) is reliable for the 6-quark, however, in employing the same 
reduction for the c-quark, one may expect quantitatively important corrections. 
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The light quark propagators are not limited in this way. They retain their full mo- 
mentum dependence, which is efficaciously characterised in the parametrisation flOfl 



a f s (x) = 2m f F(2(x + rh))) + Ffrx)?^) (b f + b f 2 F(ex)) , (7) 
2{x + m 2 ) - 1 + e- 2{x+ ™/ } 



<W) = ~ o , > ( 8 ) 

2(x + m z , u 



where: / = u,s (isospin symmetry is assumed); J-{y) := (1 — e y )/y; x = p 2 /(2D); 
fhf — rrif/y/2D; and 

4(x) := V2D4(p 2 ) , 4(x) := 2D a f v {p 2 ) , (9) 

with D a mass scale. This algebraic form combines the effects of confinement and dy- 
namical chiral symmetry breaking with free-particle (asymptotically-free) behaviour 
at large, spacelike-p 2 . The parameters: fhf, % 3 . in Eqs. (0) and @ take the values 



m 



f 



b f b{ b{ b{ 



u : 0.00897 0.131 2.90 0.603 0.185 , (10) 
s : 0.224 0.105 2.90 0.740 0.185 



which were determined flUj in a least-squares fit to a range of light-hadron observ- 
ables. The values of b[ 3 are underlined to indicate that the constraints b\ 3 = 6" 3 
were imposed. The scale parameter D = 0.160 GeV 2 . 

The heavy-quark expansion introduced above can be employed in the analysis 
of semileptonic pseudoscalar — > pseudoscalar decays: PhApi) ~* Ph 2 (P2) £ v , where 
Ph 1 represents either a B or D meson with momentum pi (p 2 = —m 2 Hl ) and P# 2 can 
be a D, K or it meson with momentum p 2 (pf> = — ^# 2 )- (Light — > light transitions 
are discussed in Ref. The invariant amplitude describing the decay is 

A(P Hl -> P H ,/v) = ^ V qQ l lfi {\ - 75 )uM^ PH2 ( Pl ,p 2 ) , (11) 

where Gf is the Fermi weak-decay constant, V q Q is the appropriate element of the 
Cabibbo-Kobayashi-Maskawa matrix (q denotes a light-quark and Q a heavy-quark) 
and the hadronic current is 

M^ PH *{ Pl ,p 2 ) := (PhMIZJuQIPhM) = U(t)( Pl +p 2 ), + (12) 

with t := — q 2 := — (pi — P2) 2 ■ The form factors, f±{t), contain all the information 
about strong interaction effects in these processes and their accurate estimation is 
essential to the extraction of V q Q from a measurement of a semileptonic decay rate. 
In impulse approximation 

M PH ^( Pl ,p 2 )= (13) 
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<Pk tr 



T H2 (k; -p 2 )S q {k + p 2 )i^S Q (k + pi)T Hl (k;pi)Sqi(k) 
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Hitherto unspecified in Eq. (|T3|) are Tn^k^pi) and r# 2 (/c;p2), the meson Bethe- 
Salpeter amplitudes. They can be obtained by solving the Bethe-Salpeter equation 
in a truncation consistent with that employed in the quark DSE. However, since 
for light-quarks we have parametrised that solution, we follow Ref. 0] and do 
the same for the light-meson amplitude; i.e., for the tt- and A'-mesons we assume 
F n ,K(k; P) = 275 (& 2 ) and employ the algebraic parametrisation [|T~0|] : 
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(14) 



m } =0 



which in concert with Eqs. (^) and (H) provides an efficacious algebraic representa- 
tion of Xn, K (k; P) := S(q + P/2) V^ K {k; P) S(q - P/2). C = 0.214 GeV is chosen 
to yield a calculated value f n = 0.131. fx = 0.160 GeV. 

For a heavy-meson, Bethe-Salpeter equation studies [O suggest the Ansatz 



^H f {k;pi) 



-)-, (l + ^iyv 



1 



<p(k 2 



(15) 



where, using Eq. 

Kf 2 - 



D, the canonical normalisation condition is 

1 iV„ 



rriHf 327r 2 



I duip(z) 2 (cts(z) + \Zua(r(zj^ :-- 



m Hf K 2 ' 



(16) 



with z — u — lE^fu and / labelling the light-quark flavour. In a solution of the 
Bethe-Salpeter equation the form of <p{k 2 ) is completely determined. However, here 
it characterises our Ansatz and we choose 



exp 



(-A: 2 /A ; 



(17) 



where A is a free parameter. As long as <p(k 2 ) is a non-negative, non- increasing 
function of k 2 , calculated results are insensitive to its detailed form. The leptonic 
decay constant in the heavy-quark limit is straightforward to determine once the 
Bethe-Salpeter amplitude is known: 



fa* 
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y / 7Jl Hf 87T 2 

from which it is clear that 



du {^/u — E) <f(z) 



<4( z ) + \\fu>av{z) 



From Eqs. (§), QT^), 



fH f V m H f = const. (19) 
and (|T9"[), and the pseudoscalar meson mass formula Hj: 
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(20) 
(21) 
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where Mu) = diag(m£, mjj, . . .) and T H is a flavour matrix identifying the chan- 
nel under consideration, it also follows |§ that 

m Hf oc m Q (22) 

in the heavy-quark limit, where rfiQ is the renormalisation point invariant current- 
quark mass. The linear trajectory becomes apparent for m# > ttik H 0- in 
contrast, for small current-quark masses, Eq. ( p0|) yields what is commonly known 
as the Gell-Mann-Oakes-Renner relation. Thus, in Eq. ( p0[ ) one has a single, exact 
formula that provides a unified description of light- and heavy-meson masses. 

Using Eqs. (|) and (pj) one finds pf from Eqs. (0) and © that the B f -> L> f 



decay is particularly simple to study in the heavy-quark limit. It is described by 
one form factor: 

1 m D ± m B 

f±(t) = ~ 2 —== £/H> (23) 



with = 1 + 2r(l - r)(w - 1), z w = u - 2EJu/W and 



(24) 



w = — J- L = v B ■ v D . (25) 

2m Bf m Df 

The minimum physical value of w is w m i n = 1, which corresponds to maximum 
momentum transfer with the final state meson at rest; the maximum value is w max ~ 
(m Bf + m 2 Df )/(2m Bf m B . f ) = 1.6, which corresponds to maximum recoil of the final 
state meson with the charged lepton at rest. The canonical normalisation of the 
Bethe-Salpeter amplitude, Eq. fllED, automatically ensures that 

£/(u> = l) = l. (26) 

Equation ([23) illustrates a general result: in the heavy-quark limit, the semilep- 
tonic decays of heavy mesons are described by a single, universal function: £/(w). 

The analysis of heavy — ► light decays is more difficult because, as remarked 
above, the current-quark mass of the u- and s-quarks m u / s < M x ~ 0(A QCD ). Hence 
the momentum-dependent modification of the dressed-quark propagator cannot be 
ignored, and the description of these decays requires a good understanding of light- 
quark propagation characteristics and the internal structure of light-mesons. The 
form factor that determines the width is 

(t) = K q> ^^- 2 F ql (t; E,m Hl , m H2 ) , (27) 

where 

4 d'y f°° 

F q >(t;E,m Hl ,m H2 ) = - / / dv / u 2 dinp{z x )£{zx) W q >(j, u,u) , 

7T J-l \/l — I 1 JO JO 



[28) 
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DATA/ESTIMATES 


/b = 0.170 GeV 


(E, A) (GeV) 




(0.442,1.408) 


V~l9 /AT 

Ir/iV 




0.48 




0.82 ±0.17 PH 


0.84 t 


/£ "(17.9 GeV 2 ) 


1.19 ±0.28 [0 


x 

1.02* 


/£* (20.9 GeV") 


1.89 ±0.53 [JlJ 


1.30* 


Br(B° -> tt^+i/) 


[1.8 ± 0.4 ± 0.3 ± 0.2] x 10- 4 |Tg 


2.0 xlO" 4 * 


— 7 \ 1 


0.18 -> 0.49 [|T§] 


0.46 


/r(o) 


0.74 ± 0.03 [17] 


0.62 
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0.91 ±0.15 ±0.06 ™ 
1.53 ±0.36 ±0.14 "-^ 


0.87 


/s. (GeV) 


0.195 ±0.035 @ 


0.184 


Ib s 1 Jb 


1.14 ±0.08 [g 


1.083 


Id (GeV) 


0.200 ±0.030 H 


0.285 


(GeV) 


0.220 ±0.030 || 


0.304 




1.10 ±0.06 [I§ 


1.066 



Table 1: Calculated results cf. data (experimental or lattice simulations) when 



we require Jb = 0.170 GeV, which is the central value estimated in Ref. [{19 
Quantities marked by * are used to constrain the parameters (E,A) by minimis- 
ing E 2 := J2iLi {[vT Xc ~ % data ]/ cr (?/)i lata ) 2 ! where N is the number of data items used. 
NB: 1) the values of jo and fjj s are obtained via Eq. ([Tj^) from f B and fs a1 re- 
spectively, using m B = 5.27, rriB a = 5.375, m D = 1.87 and m Ds = 1.97 GeV; 2) 
the experimental determination of p 2 is sensitive to the form of the fitting function, 
e.g., see Ref. |]n| and Fig. |2|; 3) an analysis of four experimental measurements of 
D s -> fiv decays yields f Ds = 0.241 ± 0.21 ± 0.30 GeV EO . 



z\ = u 2 — 2uuE, with Wq^j, u,u) depending on the light-quark propagator and its 
derivatives || . 

All that is necessary for the calculation of the mesonic semileptonic heavy — ► 
heavy and heavy — > light transition form factors, and heavy-meson leptonic decay 
constants is now specified. There are two free parameters: the binding energy, 
E, introduced after Eq. (^) and the width, A, of the heavy meson Bethe-Salpeter 
amplitude, introduced in Eq. (fiT|). The dressed light-quark propagators and light- 
meson Bethe-Salpeter amplitudes were completely fixed in the application of this 
framework to the study of 7T- and .fT-meson properties. The primary goal of this 
study is to determine whether, with these two heavy-quark parameters, a description 
and correlation of existing heavy-meson data is possible using the DSE framework. 
Some key results are presented in Table [I], which also describes how the parameters 
(E, A) were fixed. 
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Figure 2: Calculated form of £(w) cf. recent experimental analyses. The solid line 
was obtained assuming only that the 6-quark is heavy, the dash-dot line assumed the 
same of the c-quark || . Experiment: data points - Ref. f2l|; short- dashed line - linear 
fit from Ref. |8|, £(w) = 1 - p 2 (w - 1), p 2 = 0.91 ± 0.15 ± 0.16 ; long-dashed line - 
nonlinear fit from Ref. 0, £(w) = [2/(w+l)) exp [(1 - 2p 2 ) (w - l)/(w + 1)] , p 2 = 
1.53 ± 0.36 ± 0.14 . The two light, dotted lines are this nonlinear fit evaluated with 
the extreme values of p 2 : upper line, p 2 = 1.17 and lower line, p 2 = 1.89. 



The calculated form of £(w) is depicted in Fig. ||. It yields a value of p' 
—£'(w = 1) = 0.87 — 0.92,[| close to that obtained with a linear fitting form 
however, £(w) has significant curvature and deviates quickly from that fit. The 



curvature is, in fact, very well matched to that of the nonlinear fit [15 1, however, 
the value of p 2 reported in that case is very different from the calculated value. 
The derivation of the formula for £(w) assumes that the heavy-quark limit, Eq. (Bp, 
is valid not only for the 6-quark but also for the c-quark. Therefore these results 
suggest that the latter assumption is only accurate to approximately 20%; i.e., 
l/M c -corrections are quantitatively important. 

f+*(t) is depicted in in Fig. A good interpolation of the result is provided by 

= 1 °^ 58 2 > m -n = 5.67 GeV. (29) 



This value of m mon can be compared with that obtained in a fit to lattice data: f!4 
m mon = 5.6 ± 0.3. 

The calculated form of f+ K (t) is depicted in Fig. [|. The t-dependence is also 
well-approximated by a monopole fit. The calculated value of f+ K (0) = 0.62 is 



1 In this framework the minimum possible value for p 2 is 1/3 13 
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Figure 3: Calculated form of f+ n (t). The solid line was obtained assuming only 
that the 6-quark is heavy, the dashed line assumed the same of the c-quark || . The 
data were obtained in lattice simulation |14]] and the light, short-dashed line is a 
vector dominance, monopole model: /+(£) = 0.46/(1 — t/m 2 B ,), m#. = 5.325 GeV. 
The light, dotted line is the phase space factor \f+ n (0)\ 2 [(t+ - t)(t- - t)f /2 /(itm B f 
that appears in the expression for the width, which illustrates that the B — ► nev 
branching ratio is determined primarily by the small-t behaviour f? n (t). 



approximately 15% less than the experimental value [17]. That is also a gauge of 



the size of l/M c -corrections, which are expected to reduce the value of the D and 
D s leptonic decay constants calculated in the heavy-quark limit: fp = 0.285 GeV, 
f Ds = 0.298 GeV. A 15% reduction yields f D = 0.24 GeV and f Ds = 0.26 GeV, values 



which are consistent with lattice estimates [[19] and the latter with experiment f20 |. 

It must be noted that Ref. || explicitly did not assume vector meson dominance. 
The calculated results reflect only the importance and influence of the dressed-quark 
and -gluon substructure of the heavy mesons. That substructure is manifest in the 
dressed propagators and bound state amplitudes, which fully determine the value of 
every calculated quantity. That simple-pole Ansdtze provide efficacious interpola- 
tions of the calculated results on the accessible kinematic domain is not surprising, 
given that the form factor must rise slowly away from its value at t = and the heavy 
meson mass provides a dominant intrinsic scale, which is only modified slightly by 
the scale in the light-quark propagators and meson bound state amplitudes. 

This presentation illustrates the phenomenological application of a heavy-quark 
limit of the DSEs that is based on the result that the mass function of heavy-quarks 
evolves slowly with momentum. Heavy-mesons are seen to be little different from 
light-mesons: they are bound states of finite extent with dressed-quark constituents. 
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Figure 4: Calculated form of f+ K (t): the solid line was obtained assuming only 
that the 6-quark is heavy, the dashed line assumed the same of the c-quark || . 
The light, short-dashed line is a vector dominance, monopole model: f+{q 2 ) = 
0.74/(1 - q 2 /m 2 D ,), m D * = 2.11 GeV. 

The results summarised here indicate that the heavy-quark limit can be used to 
develop a quantitatively reliable description of 5-meson observables. However, it is 
inadequate for D-meson observables, where corrections of 15-20% can be expected. 
A significant feature of the DSE approach is that it provides a single framework 
for the correlation of heavy — > heavy and heavy — > light transitions and for their 
correlation with light meson observables, which are dominated by effects such as 
confinement and DCSB. 

This work was supported in part by the Russian Fund for Fundamental Research, 
under contract 96-02- 17435-a, and the US Department of Energy, Nuclear Physics 
Division, under contract no. W-31-109-ENG-38. 
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